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<N ■ Abstract 
Q ' It was recently pointed out that the fundamental Planck mass could be close to the 



0\ 



TeV scale with the observed weakness of gravity at long distances being due the ex- 
istence of new sub-millimeter spatial dimensions. In this picture the standard model 
fields are localized to a (3 -|- l)-dimensional wall or "3-brane". We show that in such 
theories there exist attractive models of inflation that occur while the size of the new 
■A • dimensions are still small. We show that it is easy to produce the required number 
of efoldings, and further that the density perturbations 5p/ p as measured by COBE 
^ ' can be easily reproduced, both in overall magnitude and in their approximately scale- 
\ invariant spectrum. In the minimal approach, the inflaton fleld is just the moduli 
describing the size of the internal dimensions, the role of the inflationary potential 
being played by the stabilizing potential of the internal space. We show that under 
quite general conditions, the inflationary era is followed by an epoch of contraction of 
our world on the brane, while the internal dimensions slowly expand to their stabiliza- 
tion radius. We flnd a set of exact solutions which describe this behavior, generalizing 
the well-known Kasner solutions. During this phase, the production of bulk gravitons 
remains suppressed. The period of contraction is terminated by the blue-shifting of 
Hawking radiation left on our wall at the end of the inflationary de Sitter phase. The 
temperature to which this is reheated is consistent with the normalcy bounds. We give 
a precise deflnition of the radion moduli problem. 



1 Introduction 



It was recently pointed out that the fundamental Planck mass could be close to the 
TeV scale 0, ^, thus providing a novel solution to the hierarchy problem for the 
standard model. Gravity becomes comparable in strength to the other interactions at 
this scale, and the observed weakness of gravity at long distances is then explained by 
the presence of n new "large" spatial dimensions. Gauss' Law relates the Planck scales 
of the (4 + n) dimensional theory, M^,, and the long-distance 4-dimensional theory, Mpi, 

Ml = (6o)"Mr2 (1) 

where bo is the (present, stabilized) size of the extra dimensions. If we put ~ 1 TeV 
then 

bo ~ 10-^^+Vcm (2) 

For n = 1, 6o ~ 10^^ cm, so this case is excluded since it would modify Newtonian 
gravitation at solar-system distances. Already for n = 2, however, 6o ~ 1 mm, which 
happens to be the distance where our present experimental knowledge of gravitational 
strength forces ends. For larger n, 1/bo slowly approaches the fundamental Planck 
scale M^,. 

While the gravitational force has not been measured beneath a millimeter, the success 
of the SM up to ~ 100 GeV implies that the SM fields can not feel these extra large 
dimensions; that is, they must be stuck on a 3-dimensional wall, or "3-brane" , in the 
higher dimensional space. Thus, in this framework the universe is (4 + n)-dimensional 
with fundamental Planck scale near the weak scale, with n > 2 new sub-mm sized 
dimensions where gravity, and perhaps other fields, can freely propagate, but where 
the SM particles are localized on a 3-brane in the higher- dimensional space. The most 
attractive possibility for localizing the SM fields to the brane is to employ the D-branes 
that naturally occur in type I or type II string theory Gauge and other degrees 

of freedom are naturally confined to such D-branes and furthermore this approach 
has the obvious advantage of being formulated within a consistent theory of gravity. 
However, from a practical point of view, the most important question is whether this 
framework is experimentally excluded. This was the subject of where laboratory, 
astrophysical, and cosmological constraints were studied and found not to exclude these 
ideas. 

There are a number of other papers discussing related suggestions. Refs. examine 
the idea of lowering the gauge- coupling unification scale by utilizing higher dimen- 
sions. Further papers concern themselves with the construction of string models with 
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extra dimensions larger than the string scale 0, ||, 0, gauge coupling unification in 
higher dimensions without lowering the unification scale [1^, the effective theory of 
the low energy degrees of freedom in realizations of our world as a brane |]1T| and radius 
stabilization |T^, There have also been many recent papers discussing various theo- 
retical and phenomenological aspects of this scheme |]13[, and a few papers on aspects 
of the early universe cosmology [III] [|T3 ITBI [ITH] that discuss issues related to those 
considered here.Q 

In this paper we will discuss infiation and other general aspects of early universe 
cosmology in world-as-a-brane models. In particular we will be concerned with aspects 
of early universe cosmology that involve the dynamics of the internal dimensions in a 
central way. We find that there exist attractive models of inflation that occur while 
the internal dimensions are still small, far away from their flnal stabilized value given 
by the Gauss' law constraint. We show that in such models it is very easy to produce 
the required number of efoldings of inflation even though no energy densities exceed 
the fundamental Planck scale, p < M^, and that supersymmetry, if it exists at all, is 
very badly broken on our brane. Further, we demonstrate that the density perturba- 
tions 6p/ p a.s measured by COBE and the other microwave background and large-scale 
structure experiments can be easily reproduced in such models, both in their magnitude 
and in their approximately scale-invariant spectrum. In the most minimal approach, 
the inflaton fleld may be just the moduli describing the size of the internal dimensions 
itself (the radion fleld of [Q), the role of the inflationary potential being played by 
the stabilizing potential of the internal space. In the case of a wall-localized inflaton, 
the cosmological constant might even result from the electroweak phase transition, in 
which case the inflaton is the Higgs. Actually, an important remark in this regard is 
that when the internal dimensions are small, b ~ M~^, the distinction between on-the- 
wall and off-the-wall physics is not meaningful: e.g., the inflationary features in V{b) 
at small b could be due to Higgs physics on the wall. 

The approximately scale-invariant nature of the primordial perturbation spectrum 



* However, we note that a closer scrutiny of ||i8|] r eveals that it is impossible to get any inflation 
at all in the models considered. The potentials in |f^ all have exponential dependence on the radion 
field, and thus lead to the power-law dependence of the scale factor on time, a ^ t^^^ |2^. The 
correct formula for the parameter (3 of is /3 = 2^2nl(n + 2), and thus j3 >2 for all cases n > 2, 
which gives only subluminal expansion, and not inflation. In fact, a complete analysis must employ 
the physically meaningful measure of scales, namely the scale factor expressed in units of the Compton 
wavelength of wall particles. (Because of the possibility of Weyl rescaling between string and Einstein 
frames one must be careful about the definition of physical quantities.) In such physical units it is 
possible to show that the scale factor always expands sub-luminally, even more slowly than the naive 
Einstein frame expansion, and thus none of the solutions of llq] contain a stage of inflation. 
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implies that, during inflation, the internal dimensions must expand more slowly than 
the universe on the wall. Thus we are led to consider a form of asymmetric inflationary 
expansion of the higher-dimensional world. It is very interesting to note that requiring 
the consistency and naturalness of both the long duration of inflation {i.e., number of 
efolds Ne ^ 100), and the magnitude of Sp/p, implies that inflation should occur early, 
with the internal dimensions close to their natural initial size bi ~ M~^. In particular, 
we avoid the introduction of extremely light inflatons which seem to be needed if 



inflation occurs after the internal dimensions reach their current, large, size |T^, |T5|. 
We emphasize here that this constraint emerged from the assumption that during 
inflation the internal space was already large and stable. In fact, in such scenarios 



obtaining 6p/p requires that the inflaton is even lighter [|T6|, and further, that inflation 
occurring only after stabilization cannot explain the age of the universe |16]. Indeed, 
the wall-only inflation cannot begin before t ~ H^^ ~ Mpi/M^ » M~^, when the 
universe is already very large and old. 

However, in the models of inflation at very early times which we consider, the effective 
4D Planck mass is both much smaller than it is now, and in general, time-dependent 
due to the variation in the volume of the internal dimensions. Moreover, inflation 
occurs soon after the "birth" of the universe when the sizes of all dimensions are close 
to their natural initial size ~ M~^. Hence early asymmetric inflation solves the age 
problem too. 

The framework in which we discuss these issues is that of semiclassical (4 -|- n)- 
dimensional gravity with an additional potential V{b) depending on the size of the 
internal dimensions. As discussed in Ref . this rapidly becomes a good approximation 
at energy scales below M^,, which, self consistently, is the correct regime for early 
inflation, essentially because of the COBE constraint that the density perturbations are 
small 6p/p ~ 2 X 10~^. Most of our analysis will utilize the so-called string frame with 
explicit scale factors (a(t), b{t)) for the sizes of our brane- localized dimensions and the 
internal dimensions respectively. In this frame the effective long-distance 4-dimensional 
value of Newton's constant depends on b{t) and changes with time. However this is 
the frame where the "measuring sticks" of particle masses and Compton wavelengths 
are fixed. On the other hand, in the analysis that follows we will find it occasionally 
useful to employ the correspondence of the low-energy (4 -|- n)-dimensional theory with 
an effective 4-dimensional theory which takes the form of a scalar-tensor theory of 
gravity, the radius of the large internal space b playing a role similar to the Brans- 
Dicke scalar field. The inflationary dynamics from the point of view of this frame 
is completely equivalent to the usual slow-roll scenarios, and the conditions for the 
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asymmetric expansion of the universe on the brane relative to the internal dimensions 
are completely equivalent to the slow roll conditions in the usual inflationary models. 
Further, this analogy provides very useful consistency check for the determination 
of the density perturbations. The relationship between the scale factor b(t) and the 
canonically normalized 4-dimensional ("radion") fleld is given in equation (^). (We 
will often abuse our terminology and refer to b{t) itself as the "radion".) 

It is important to realize that under quite general conditions, the early inflationary 
era is followed by a long epoch where the scale factor of our brane-universe undergoes 
a slow contraction while the internal dimensions continue to expand towards their flnal 
stabilized value. We show that even with the inclusion of a potential for b, it is possible 
to exactly describe the evolution during this epoch, and we present a class of exact 
solutions which generalize the usual vacuum Kasner solutions. The total amount of 
contraction of our universe on the brane is bounded by a small power of the expansion 
of the size of internal dimensions, and varies between at most 7 efoldings in the case of 
two extra dimensions to at most 12 efoldings when there are six extra dimensions. We 
show that during this phase of b{t) evolution to the stabilization point, the production 
of bulk gravitons by the time- varying metric remains completely suppressed, ensuring 
that the bulk is very cold at, and after, the stabilization of the internal dimensions. 

However, the particles produced by gravitational (Hawking) effects on the brane at 
the end of the inflationary de Sitter phase can play an important role. Their energy 
density is blue-shifted by the slow contraction of a{t), and when this wall- localized 
energy density exceeds the energy density in the radion, the contraction of a{t) ceases, 
with a "Big Bounce" occurring, leading on to an expansion of both a and b. It is very 
tempting to view the reheating on the wall as due to this contraction of primordial 
de Sitter radiation. This is especially so as the flnal temperatures come naturally close 
to the normalcy bounds (the upper bound T* on the temperature of 0] that ensure 
that the evaporation to bulk gravitons remains negligible). 

However, such a conclusion is premature, since as we argue below the energy density 
in the radion is naturally of the same order as that on the wall, and this leads to a 
radion moduli problem in the later evolution of the universe. Explicitly, because of 
the gravitationally suppressed couplings of this fleld to wall-localized SM states, the 
radion decays very slowly back into light wall flelds. As a result, its energy density 
would red-shift only according to pb ~ l/a^ compared to p^uaii ~ and comes to 

dominate the total energy density. This is just the standard moduli problem, which is 
inextricably linked with the mechanism of reheating in world-as-a-brane models. 
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Concretely, the basic picture we advocate is: 



The quantum creation of the universe takes place with the initial size of all 
dimensions close to the fundamental Planck scale M~^. In particular the initial 
size bi of the internal dimensions is of this size. 



• A prolonged period of inflation in a direction parallel to our brane takes place, 
with our scale factor a{t) increasing superluminally a{t) ~ with p » 1, and 
b{t) — hi essentially static. (In the cases that we explicitly discuss the inflation is 
quasi-exponential, a{t) ~ exp(ift).) This brane inflation is driven by either the 
stabilizing potential of the radion b itself, or by a wall-localized field with effective 
non-zero cosmological constant. In either case it is unnatural for the size of the 
effective 4-dimensional cosmological constant to exceed roughly ( TeV)^. Since 
the internal dimensions are small the effective 4-dimensional Newton's constant 
is large 

Gn, initial = " J^' 

Thus the Hubble constant during this initial period of inflation can be large even 
though the energy density is quite small, {V) ~ O(TeV^), 

(4) 

mfl l,n^n+2 V ) 

Thus inflation can be rapid, and moreover, as we will argue in detail the density 
perturbations can be large, being determined in order of magnitude to be 

p - M,{M,biYl^S' ^ ' 

where 5" is a parameter that encapsulates both the duration of a{t) inflation and 
the deviation of the perturbation spectrum from the scale-invariant Harrison- 
Zeldovich spectrum. (We will argue that S < 1/50.) We can turn this around by 
imposing the COBE-derived normahzation on 5p/ p, and thus on (V) /6"M"+^. 
During this period the size bi of the internal radii are flxed by the over-damping 
arising from H^n. 

• Wall inflation now ends, with H starting to drop and simultaneously the radion 
starting to evolve to its' minimum at 6o- The initial motion of the radion can be 
quite involved, but from the coupled equations of motion for the scale factors a{t) 
and b{t) we will see that as b evolves towards its' minimum, our scale factor, under 
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very general conditions, undergoes a collapse. The detailed behavior generalizes 
the well-known Kasner solutions where some dimensions expand, while others 
collapse, both with determined (subluminal) power-law dependence on time. 

• The contraction epoch ends when the blue-shifting radiation density on our brane 
becomes equal to the energy density in the radion field pwaii = Pb- A model- 
independent initial source for the radiation on the wall is the Hawking radiation 
left over from the early infiationary de Sitter phase. After the contraction of a{t) 
reverses, the radion and wall-localized energy densities scale together until the 
stabilization point is reached. 

• Finally around the stabilization point 60 the radion field starts to oscillate freely. 
Since this energy density scales as l/o?, and the wall-to-radion energy densities 
are initially comparable at the start of the oscillation era, the radion energy starts 
to dominate the total energy density. 

• The most serious question that early universe cosmology presents in the world- 
as-a-brane scenario is how do we dilute this energy in radion oscillations to an 
acceptable level. The radion is long-lived, its' decay width back to light wall 
states being given byQ 

r, ^ (6) 

We thus require some dilution in the radion energy density, either by a short 
period of late inflation followed by reheating, or by a delayed reheating after pb 
has sufficiently red-shifted. The amount of dilution of the radion energy density 
that we require is given roughly by leV/T^, ~ 10~^, so that only about 5 efolds 
of late inflation is needed. 

The paper is organized as follows. In Section 2 we will discuss the details of the 
inflationary stage, both in the original string (brane) frame and in the effective Einstein 
frame, showing the initial conditions, the conditions on inflationary potentials and the 
resolution of the COBE constraints. Section 3 is devoted to the study of the era of 
radion evolution to stabilization. In subsection 3.1 we derive the exact solutions which 
generalize simple vacuum Kasner solutions, and which demonstrate that a period of 



^Here m^p is the mass of the canonicaUy normaUzed field corresponding to h. The experimental 
bound on this mass is rriip > 10~^ eV. Note that the decay width is increased if we take into account 
the possible presence of many branes in the bulk. Indeed some scenarios of radius stabilization e.g., 
the "brane crystallization" picture [Q, require -/Vwaii — (Afpi/M*)^^"~^)/" branes in the bulk. If each 
of these have 0(1) light modes then the total decay width to all branes is greatly enhanced. 
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a(t) contraction subsequent to inflation is quite generic. In subsection 3.2 we discuss 
the phenomenological consequences of tliis era, including the cause of its cessation and 
the (minimal) production of bulk gravitons. There we will also show how the normalcy 
bounds for the temperature of our brane after stabilization are met by the evolution. 
In subsection 3.3 we discuss we discuss the issue of the production of bulk gravitons 
during the era of radion evolution and show that the quantities produced are harmless. 
In subsection 3.4 we provide a physically quite useful discussion of the contraction and 
stabilization eras in the Einstein frame, where many of the arguments of the previous 
subsections can be understood quite simply. In subsection 3.5 we provide a short 
discussion of the radion moduli problem, and quantify its' size using the evolution 
equations discussed earlier in Section 3. Section 4 gives our conclusions. We also 
include three extensive Appendices: Appendix A contains a short discussion of some of 
the basic kinematics of brane evolution embedded in higher dimensions. Appendix B 
gives some additional details of the exact solutions for the post-inflation slow evolution 
era where a{t) contracts. Appendix C proves that the contraction of a{t) is self-ceasing 
by the blue-shifting of brane-localized radiation, and also discusses the exact solutions 
that may be obtained for the "Big Bounce" that terminates this contraction. 



We now embark on our detailed discussion of the evolution of the sizes of our brane 
and the transverse "internal" dimensions. In particular in this section we focus on the 
physics of the (early) inflationary epoch. 

We are interested in the case of a 3-brane embedded in a (4 + n)-dimensional space- 
time. To reproduce our observed world, ultimately the 3 spatial dimensions parallel to 
the brane must be as least as large as our current horizon size, while the n transverse 
spatial dimensions have to stabilize at a size bo given by the constraint on Gn, (|I|), 
(6o)"M"+^ = Mpj. Note that we take the internal space to be topologically compact 
from the beginning. That is we impose periodic boundary conditions in the directions 
transverse to the wall. These conditions reflect the fact that the low energy theory is 
four- dimensional, and that the classical evolution preserves the topology. 

The total action is comprised of a bulk part. 



2 Early inflation 




(4+n) 




■matter 



) 



(7) 



7 



and a brane part, 

5'brane — J d X^J det Q^^^^'^'^^ ^-^standard model 

where ^matter is the Lagrangian of the bulk fields apart from the graviton. These fields 
give rise to the stabilizing potential V that we discuss below. The ellipses denote 
higher- derivative terms that can be safely ignored in most of the regime of interest 
as curvatures are small compared to the fundamental Planck scale M^, apart possibly 
from the very early pre-infiationary stage immediately after the quantum creation of 
the universe. However any signatures of such a phase of high energy and large curvature 
have been wiped out of the visible universe by the subsequent stage of infiation. We 
will therefore ignore this stage as practically invisible, and begin the description of 
the world-as-a-brane universe at curvatures an order of magnitude or so below the 
fundamental Planck scale. Indeed, at and below these scales, the description based 
on the actions (|^) and should be reliable. The background metric for the (4 + 
n) -dimensional spacetime which is consistent with the symmetries of the brane-bulk 
system is of the form 



( ^ \ 

-a{tfgij 



(9) 



where a is the scale factor of the 3-dimensional space, and b is the scale factor of the 
internal n-dimensional space, with geometry set by Qij where det((?jj) = 1. 

As shown in Appendix A, the equations of motion for the coupled {a(t), b{t)} system 
can be written in the form. 

Ha + nHaH, + ^Hl = J— (b^ - (n - 2)1^) (10) 

2(n + 2)6"Mr+^^ V J 



and 



H, + 2>HaH, + nHl = ^ ^ . , f-^ - 2v\ (11) 



[n 



where we have introduced the Hubble parameters Ha = a/a and Hb = b/b for the two 
scale factors, and an overdot denotes a derivative with respect to t. These equations 
of motion are supplemented by the constraint equation 

6Hl + GnHaHb + n{n - l)Hl = ^ (12) 



8 



Note that in these formulae the potential V is the effective 4- dimensional potential with 
mass dimension [\^] = 4 - i.e., that arising from projecting the bulk energy density 
that results from i2mattcr onto our wall. 

The natural inflationary initial conditions consist of taking the initial (4 + n)- 
dimensional universe to be relatively smooth, flat and potential energy-dominated 
over scales roughly given by the fundamental Planck length and not much 

greater In the brane language the initial conditions we adopt mean that 



we take a portion of the brane of linear size ~ embedded in a (4 + n)-dimensional 

volume of the same linear size, to be relatively flat, smooth and straight. Indeed, it is 
most natural to suppose that the universe starts as a small (4 + n)-dimensional domain 
of linear extent in all directions, and then undergoes an epoch of inflationary 

expansion almost immediately, at least in the directions parallel to our brane. 

Quantitatively these conditions can be stated in terms of the initial conditions for 
the brane and bulk horizons as follows. The size of initially isotropic and homogeneous 
causal domains must be 

H-,l ~ H^,l > M;' (13) 

Moreover we can also define the initial conditions for the scale factor of the brane a and 
the radion field h by referring to the flatness problem. Since close to the fundamental 
Planck scale the energy density on the wall and in the bulk is of order unity in funda- 
mental Planck units, if we pick the gauge such that any intrinsic spatial curvature on 
the brane or in the bulk is k = ±1, it cannot exceed in magnitude Pp, leading to the 
conditions that[] 

a. ~ &i > (14) 



It is easy to recognize the conditions ([13| ) and (|l^ as precisely the consistent initial 
conditions for inflation, in the case when the fundamental Planck scale is M^. How- 
ever these initial conditions are merely a rough estimate coming from the requirement 
that the approximation based on semiclassical gravity is valid, and that the usual cos- 
mo logical problems (horizon, flatness, homogeneity etc..) are solved by subsequent 
inflation. 

Further, the generation of sufficient density perturbations Sp/p ~ ^p/p|cobe ~ 10~^, 



^In spatially flat FRW universes the magnitude of the scale factor a is physically meaningless. 
However, the spatial curvature, namely the quantity k/a^ is physical. Thus choosing the value of a 
at some instant t corresponds to specifying curvature of spatial hypersurfaces. Equivalently, choosing 
the value of the constant fc, by normalizing it to ±1, corresponds to picking the units for the scale 
factors. We will assume this throughout this work, even if we do not explicitly specify the spatial 
curvatures (which can be ignored after prolonged inflation). 
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implies, at least without the unappetizing introduction of very light fields on our wall, 
or of the assumption of very unusual initial conditions, that the inflationary epoch that 
solves the flatness and horizon problems occurs when the size of the internal dimensions 
is still relatively small. The basic reason for this is that in a theory with fundamental 
Planck scale M^,, the energy density localized on the brane should not exceed M^. 
Similarly the energy density in the bulk should not exceed M""*"^. Moreover, the 
(matter) energy density of the universe at its' birth is expected to be of order M^"*"^. 
If inflation occurs when the internal dimensions are already large, then the effective 
Newton's constant in our 4-dimensions is already very small GN,efr = 1/{M^^'^IP') <C 
and energy densities of order or less will lead to a very small Hubble constant, 
and thus typically unacceptably small density perturbations.| 

To see this in detail consider the expression for the density perturbations generated 
during slow-roll inflation driven by a (canonically normalized) fleld 

p 12tt if ' ^ ^ 

(The use of this 4-dimensional expression will be self-consistently justifled later on in 
our analysis.) The Hubble constant on our brane during inflation is given in order of 
magnitude by H[^f{ ~ l^/(M"+^6"), where 6/ denotes the size of the internal dimensions 
during the inflationary epoch. Since V < if ^ 1 then 6p/p is very small unless 
if is extremely small. Although this is a logical possibility, it requires extraordinary 
flne-tuning and we will not consider this case here. On the other hand if is 0{1) 
then sufficiently large density perturbations easily result. Moreover, from the power 
dependence of ifinfi on bj we see that if the size of the internal dimensions changes 
significantly during inflation the spectrum of density perturbations will be very far from 
scale invariant. Since the spectral index of density perturbations Up is constrained by 
the cosmic microwave background (CMB) and large scale structure measurements to 
be not significantly different from the scale-invariant value = 1, 

|np-l|<0.2, (16) 

it is necessary that the evolution of the internal dimensions is slow compared to that 
of our scale factor a: Hb <C Ha. A successful phenomenology thus results if the ratio 
Hjj/Ha approaches a zero for a (small) range of b around the value bj. Let us expand 

§Note though that if the internal dimensions are large b ^ ^/M■^, a small bulk (4 + n)-dimensional 
energy density still allows the effective "projected" energy density on the wall to exceed M^, without 
requiring the use of the full quantum theory of gravity; the semiclassical approximation is still good. 
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the ratio around this point 

( h \^ 

(17) 



Hi 
H, 



a 




The dimensionless parameters S and T will be bounded in size by the requirements 
that the spectral index n of the density perturbations is close to 1, and that sufficient 
e-folds of inflation occur to solve the horizon and flatness problems. Furthermore Hi, 
must be small as compared to H^, otherwise the spectral index would again too quickly 
deviate from 1. 

To quantify these restrictions let us return to the equations of motion for a{t) and h{t) 
and the general expression for the density perturbations. Successful inflation requires 
that the number of efolds of superluminal expansion be sufficiently large. If the 
initial and final values of the scale factors just before and just after the inflationary 
epoch are denoted by subscripts (i, /) respectively, then 

iv. = r"^ (18) 

h Ha dh 

h 1 dh 



S + T{h/hi - 1)2 h 

To a good approximation inflation begins and ends when the condition Hb/Ha ^ 1 is 
violated. Hence to obtain inflation, there are two possibilities. Either S must satisfy 
<^ 1, which in turn implies 

hl = l^T-'/\ (19) 
hi 

The alternative is to have 5* small (but not excessively) while T must be extremely 
small, which corresponds to slow, power law, inflation. Substituting these endpoint 
values into (^) and performing the integral gives the constraint on S and T arising 
from Ng. 

Requiring, say, > 100 then puts an upper bound on the size of (5*, T). In various 
limits this relation becomes easy to state explicitly. For example if S" = T -C 1 then 

Ne^^ (21) 
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while if S is very small, while T remains 0(1), 

iV,^-^-log(l/5). (22) 

In any case obtaining ~ 100 efolds of inflation requires either S" ~ T < 0.02, or 
T ~ 0(1), and 5* < 10"'^. So in practice to get a sufficient number of efolds S should 
be somewhere in the range 0.001 < S < 0.02 (or, of course less) depending on T. 

Now let us turn to the constraint arising from the magnitude and spectral index 
of the density perturbations. To do this in detail, we must make some assumptions 
about the identity of the infiaton. Since we have argued that inflation should occur 
far away from the eventual stabilization value of b, the most natural candidate for the 
potential energy that drives inflation is the radion potential itself, with therefore the 
radion playing the role of the infiaton. The correctly normalized field ip corresponding 
to the scale factor b is given by (see Appendix A) 



= ,j2n{n - l)Mi"+2)/25(n-2)/2^^ ^23) 

Substituting this into the expression for the density perturbations ([l5|) we find 

6p 5 H', 



P ~ 12n^2n{n - 1) Mi"+')/'6("-2)/26 

5 Hg 

l2^j2n{n-l)M.{M^hYl^S' 



(24) 



where in the last step we have used (|T^) and the fact that we are close to 6/. Requiring 
b{t) to be essentially static during inflation, and in particular that Hi,/ Ha and H^/H^ 
are small near bi is equivalent, using the equations of motion, to the statement that 



2 (^dbV - 2V 



3 {bdbV - (n - 2)V) 

which translates into the condition that 

fb 



< 1. (25) 

b~br 



-dbV - 2V 
n 



« V{bj), bjd,V{bj). (26) 

bi 



Therefore {bdbV)\bj — 2nV\bj, and up to small corrections, the Hubble parameter Ha 
during inflation is given by 

H^ ~ (27) 
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Finally this together with (p5[) gives the expression for the inflation generated density 
perturbations 

Sp 5 1 f V{bi) 



1/2 



p Un S{M^hiY \l2n{n - l)Mt J 
The spectral index Up is defined by the comoving wavenumber dependence 

6p 



^K-i)/2 



P 



(28) 



(29) 



where at horizon crossing we have the relation k/a = Ha- From (|29| ) we can extract a 
more convenient expression for Up-. 



rip -I 



^ d\og{5p/p) 
d\og{a) 
d\og{6p/p) 
d\og{b) 



2S 



where in the second line the parameterization for Hb/Ha ([T7|) has been used. Applying 
this formula to the expression for the density perturbations ( pSD leads to 



n. 



l = 2S\-n + 



dV 



2Vd\og{b) 



(31) 



But from the two expressions for the ratio Hh/Ha, one from the equations of motion 
(|iy,|l2), and the other given by the parameterization (|1^ we have, after some algebra, 



dV 



and therefore 



2Vdlog{b) 



n-^^^{S + T{b/b, 



n{n + 2) 



(S^ + ST{b/bi 



(32) 



(33) 



This expression together with the experimental constraint |np — 1| < 0.2 forces S < 0.1, 
thus ruling out the case of T <^ 1, S' ~ 0(1) allowed by the earlier condition of a 
sufficient number of efolds of inflation. Note however that the solutions with T ~ 0(1) 
and S < 0.01 automatically satisfy |np — 1| < 0.2 over essentially the whole range of 
a{t) inflation. (Experimentally all we require is that |np — 1| < 0.2 in the range of 
scales between the COBE and large-scale structure measurements - roughly 10 efolds 
rather than the full duration of inflation.) 

Finally we can use the magnitude of the measured CMB fluctuations to constrain 
the size bj of the internal dimensions during the inflationary epoch. Recall that COBE 
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and the other CMB measurements tell us that at the time the scales k ~ 7iJnow were 
being inflated outside of the horizon, the density perturbations were of size 

^ ~ 2 X 10-^ (34) 

P COBE 

This together with the formula for Sp/p (pHf), and the bound on the parameter 5* leads 
to a constraint on a combination of the size of the potential during inflation V{bi)/M^ 
and the volume V} = (6/)" of the internal dimensions, both in units of the fundamental 
Planck scale: 

For 5 ~ T ~ 0.02 and V{bj) having a perfectly reasonable value of V{bj) ~ (200 GeV)"^, 
(for example if M^, ~ 1 TeV) we therefore discover that the epoch of inflation gener- 
ating COBE needs to occur around the value 6/ ~ 10'^/" M~^, i.e., when the internal 
dimensions are still relatively small as expected. 

Now let us return to the justiflcation of the use of the usual 4-dimensional expression 
(p!5|) for the density perturbations. If the de Sitter horizon {Ha)~^ during inflation is 
much smaller than the size of the internal dimensions then the full (4 + n)-dimensional 
expression for the density fluctuations must be employed, while if (Ha)'^ ^ bj then 
it is correct to use the long-distance, effectively 4-dimensional, expression. From the 



constraint arising from 6p/p (|35|), and the expression for Ha during inflation (27) we 
flnd that 

(5 X 10-3)("-2)/" / S x("-2)/2n /^^^^^X (n+2)/4n 



Thus Habj is always substantially less than 1 for the range of parameters of interest, 
and the 4-dimensional description of the generation of Sp/p is correct. 



2.1 Inflation from the Einstein frame perspective 

Most dynamical aspects the world-as-a-brane scenario are the most transparent in 
terms of the geometric variables employed so far, which can be referred to as the 
string-frame quantities. The reason is that the kinematics in this frame is automatically 
expressed in terms of the units felt by the observers which live on the wall. In particular, 
such observers choose to deflne scales using their own masses and Compton wavelengths 
as yard sticks, and in the string frame these yard sticks are time independent. However, 
it is illustrative to consider geometrical evolution of the universe in light of the reference 
frame where the gravitational sector of the theory coincides with classical Einstein 
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gravity. In this frame, the equations of motion resemble the coupled gravity-matter 
models considered in the context of non- minimal theories of gravity , and the results 
obtained from them can be a useful addition to physical intuition. 

For this reason, here we will recast the description of the inflationary stage into the 
Einstein conformal frame of the theory on the brane (|^), This frame is defined 
by the requirement that the propagator of the graviton is canonical. The map which 
casts the kinetic term of the graviton into the canonical form is 

f M V 

V = [j^J {M^hTg,^ (37) 

Explicitly, we can define the Einstein frame comoving time and scale factor according 
to 

<rt = ^{M^br^'dt a = -^{M^hY'^ (38) 

ikZpi ikZpi 

where all barred quantities refer to the Einstein frame. Now, the radion field is equiv- 
alent to a scalar field, defined by 



M-V^M5exp(-^^^] (39) 

Note that the normalization employed here is dictated by the definition of the pertur- 
bation of the radion away from its mean value during inflation, which sets the effective 
background Planck mass. We note that the initial conditions for inflation in the Ein- 
stein frame, in the units of the effective Planck scale, require homogeneity and flatness 
over distances of ~ lpi,eff = 1/^pi, which can be recognized as a usual inflationary 
initial condition. 

This Einstein frame picture is very useful to compute the density contrast in models 
where radion is the inflaton. The density contrast is given by the standard formula 

6p h'^ 

— = = 8 — 40 

in our normalizations, and where the prime is the derivative with respect to t. To 
determine the density contrast in the string frame, we recall that it is, roughly, con- 



formally invariant during inflation [22|, and conformally transform (HQ) to the string 



frame. The radion and the Hubble parameter transform to 



if' = -J2n{n - 1) 



M,(M,6)"/ 
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which is straightforward to determine from their definitions. With this, we find 

5p _ 8 {Ha + ^H,f 



(42) 



Now we consider the tilt of the perturbation spectrum. To consider it, we define the 
effective slope by 

where the last equahty arises because of the slow roll conditions Ha,b{ta) ~ Ha,b{h) 
during inflation. Since h{ta) ~ h{tb){l-\-Hi,5t) and since we compare the tilt between 50 
and 60 efoldings, we get 5t ~ 10i?~^. Thus roughly, ^ ~ 1 — On the other hand, 

/, \(np-l)/2 

during this time 9 changes according to ^ ~ ( |^ j , and given that ~ Ofei^a, we 

find that the ratio of the wave vectors is given by ka/ k^ ~ aa/ ~ exp{Ha5t) ~ exp(lO) . 
Taking this and the bounds on the spectral index Up, < 1 ± 0.2, we find that 

e ~ e±^/2. (44) 

This leads to the following inequality which the expansion rates must satisfy during 
infiation: 

1 — exp(±l/2) 
As a result, the density contrast 6p/p is 

Sp^, 8 H', 



(46) 



The COBE data tell us that 5p/ p ~ 10^^. Using this and Ha > 15nHb, we find 

Ha < 10-^M,(M,6)"/2 (47) 

We can use this inequality to obtain a bound on the size of the internal dimensions bj 
during inflation. We repeat that the evolution is similar to the Einstein frame thanks 
to the slow roll conditions. Hence the vacuum energy density during inflation is 

Pvacuum ~ M^^+^'^b^jHl (48) 

Using our estimate for Ha, 

lQ-^'Mt{M,bif" (49) 
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However Pvacuum must not exceed M^. In fact, since initially the theory is very close to 
the quantum gravity scale, in order to insure the validity of the semiclassical approxi- 
mation the energy density cannot exceed, let's say, 10~^M^, meaning that the energy 
scale is about a factor of 10^/^ below the quantum gravity scale. This ensures that the 
semiclassical description is correct. Using the equation for the energy density in terms 
of bi above, we obtain the inequality 

bi < lO^/'^M-i (50) 

Note that the upper limit on bj ranges between about 200M~^ for n = 2 to about 
6M~^ for n = 6. In either case however, these are clearly the correct initial conditions 
for inflation, which in fact come naturally in this context. 

3 Post-inflation evolution to stabilization point 

As we have discussed in the previous section, when the inflationary stage ends, the ratio 
of internal space to on-brane Hubble rates of expansion HiJ Ha approach unity. From 
this time on, the slow roll conditions for the effective potential cannot be upheld any 
more. The kinetic energy stored in the expansion of the internal space begins to play 
a significant role in the evolution of the brane-world. Rather interestingly, however, 
we will see that under rather general conditions the expansion of the whole higher- 
dimensional universe under the combined influence of the potential and kinetic energies 
during this era is well approximated by a generalization of the well-known Kasner 
solutions. These solutions generically describe cosmic evolution which is anisotropic 
in different directions, with a subset of the directions contracting while the others 
expand. Specifically we will find that after the end of the initial brane-world infiation 
the directions longitudinal to the brane contract, while the internal directions expand, 
both according to some power-law time dependence. The precise form of the power 
laws is controlled by the codimension of the brane {n in our notation) and the leading- 
order behavior of the stabilizing potential V{b) as a function of the scale factor b of 
the internal dimensions. 

3.1 Theory of the era of contraction 

To illustrate generic features of such behavior, and show why it should be immediately 
expected, at least in a subset of cases, consider the limiting case where after the exit 
from the infiationary stage the potential V{b) drops by many orders of magnitude. 
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The evolution of the scale factors a{t) and b{t) is then controlled entirely by the kinetic 
energy, and the equations of motion appropriate for this case are just the higher- 
dimensional vacuum Einstein equations, 

R^.u = (51) 
In other words ( pjOjpTD with V set to zero. The (4 + n)-dimensional solutions which are 



consistent with the brane symmetries are of the form 

ds-' = -dt^ + al(^^^ dxl + bl(^j^ dyl (52) 

where Oj, hi are now the "initial" values of the scale factors as set by the end of the infla- 
tionary stage. The powers k and C are uniquely determined by the Einstein equations, 
which reduce to two simple algebraic equations, 

3A; + n£ = 1 
2,k'^ + nf = 1 (53) 



A metric of the form (|52[) with exponents satisfying ([53| ) is known as the Kasner 
solution. The solutions of the algebraic equations for the exponents are 



3^ ^3n(n + 2) 
3(n + 3) 



n ± J3n(n + 2) 

^ = , ' , (54) 

n[n + 6) 

Phenomenologically, we certainly need the internal dimensions to grow in size and 
approach the stabilizing value, and this selects the upper sign in the equalities (0), 
which then implies that our longitudinal brane directions contract. As we will soon 
see when we re-introduce the potential V{b), this behavior is physically selected by 
the asymmetry embodied in the potential, and in particular the fact that it must have 
a minimum at the stabilizing value b^, with 6q = Mpj/M""*"^. In any case in this 
simple potential-free Kasner case, the explicit values for the powers range between 
k = -0.1266, £ = 0.69 for n = 2, and = -1/3, £ = 1/3 for n = 6. 

It is useful to express the brane scale factor a{t) as a function of the size of the 
internal dimensions b{t), 



a,; V bi / 



(55) 
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The amount of contraction of our scale factor is thus determined by the increase in b 
(which in turn is determined once the end of inflation value of 6 = 6j is specified), and 
the known power —\k\/i. Note that for n = 2, k/i ~ —0.183, and so the brane world 
contracts by only one order of magnitude for almost six orders of magnitude of radion 
increase. As n increases, this dependence speeds up, and for n = 6 the contraction of 
the brane dimensions and expansion of the internal dimensions are essentially equal in 
magnitude (we have of course ignored small transient effects at both the start and end 
of this Kasner phase which slightly modify the above relationships). 

We are typically interested, however, in the case where there exists a non-trivial V{b) 
potential as well as matter on our brane. We again take the metric to be of the form 

ds"^ = -dt^ + a^{t)dxl + b'^{t)dyi (56) 

and assume that the matter on the wall is represented by a classically conserved, perfect 
fluid energy-momentum tensor T^^ = (p + p)u^Uu + pg^u, V^T^'^ = 0, where is a 
future-oriented timelike vector, with the components = (l,03+„) in the basis (0). 
Here p is the energy density of the wall matter and p = 7p is the pressure, with 7 a 
constant given by the speed of sound on the wall. The equations of motion slightly 
generalize (|^) and (|ll]). 



6 ' ' " M:^+26"\^n + 2 n{n + 2) db 2{n + 2)J 

'l.OH^^vHH _ ^ ( b dV n-2 p + (n - l)p \ 

a ^ " 1^2(n + 2) 96 2{n + 2) 2(n + 2) ) 

p + 3Ha{p + p) = 0. (57) 

The last equation is of course the usual wall energy-momentum conservation equation, 
V^T^"" = 0. 

It is important to note that the contraction of the brane in Kasner-like solutions 
leads to an increase of the energy density of any matter or radiation present on the 
brane at the beginning of the contraction, and this will lead to "bounce" solutions 
that will be important to us later. Also note that these equations generically receive 
quantum corrections from particle production via curved space effects or the conformal 
anomaly. The corrections would manifest themselves as a nonzero source term on the 
RHS, which as we will see would lead to significant effects only early on, and therefore 
could be modeled by an appropriate choice of the initial condition for p. 
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Namely, the energy density on the wall corresponds to the particles produced by 
changing gravitational fields at the end of infiation. Since these phenomena are es- 
sentially similar to Hawking radiation, the value of p as compared to the potential is 
suppressed by a factor of H'^/M^^'^IP at the end of inflation. Using the COBE data to 
constrain this quantity, we see that initially p is much smaller than V . Hence ignoring 
it is an excellent approximation. For most of this stage, therefore, the dynamics of 
the universe is determined by the interplay of the radion kinetic and potential energy, 
which completely determine the evolution of the wall geometry. We will simply carry 
on with the analysis of (|57|) , treating p as a small perturbation and ignoring its back- 
reaction on a and b. However, although to the lowest order, the particle production 
phenomena are negligible, they could lead to interesting effects for reheating. We will 
return to this later. 

However before we consider such blue-shifting of the brane-localized energy density, 
let us turn our attention to the more general and appropriate case of non-negligible 
radion potential. Remarkably in this case a form of Kasner-like behavior still applies. 

Quite generically, in the semiclassical limit, fully valid at this stage of the evolution, 
the potential may be viewed as an expansion in inverse powers of the radion field b. At 
generic values of b away from the stabilization point the potential will be dominated by 
a single term in this expansion. Hence, we can simply approximate V{b) by a monomial 
of the form V = Wb~^, where is a dimensionful parameter with [W] = 4 — 

If we substitute V = Wb~^ in the equations of motion, ignoring p and p for the 
moment we find 



These equations can in fact be solved exactly! With appropriate field redefinitions 

and gauge (coordinate) transformations, they can be mapped to a system of equations 

describing the motion of two particles in one dimension, one free, and another Liouville 

with an exponential potential. 

First note that in a certain parameter range these exact solutions asymptotically 

could well depend logarithmically on h. This mild additional b dependence will not change 
either our qualitative, or to a good approximation our quantitative conclusions. For simplicity we 
ignore it in the following discussion. 





(58) 
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converge to the "potential-free" Kasner solutions ( p2D and (Q) above, as we discuss in 
Appendix B. This parameter range turns out to be the one in which, upon substitution 
of the power dependence ( |5^ ) into the equations ([58|), the potential terms on the RHS 
vanish more quickly as a function of time than the LHS. Since the LHS always scales 
as ~ this is the case when (n + p)^ > 2, or equivalently 

{n + p) 



„,„^3).n + \/3n(n + 2)>2 (59) 

This gives a curve in {n,p) space above which {i.e., for larger values of p) the exact 
solution asymptotes to the potential-free Kasner exponents. The critical values of p 
vary from 0.899 at = 2 to p = at n = 6. 

We now discuss the changes of variables which allow the exact solution of the equa- 
tions with potential. First define 

a = a,e"(*) b = he^^'^ (60) 

where Oj and hi are as before the initial values of the scale factors at the beginning of 
the epoch of radion coasting. This means that the appropriate initial conditions for 
the dimensionless variables a and (3 are ai = f3i = 0. Further, defining the parameter 
uj = n+2. n+p , substituting this and ( PH ) into (|^), and going to the new time variable 

* i 

t: 

dr = -e-^'^-'^'^dt (61) 

leads after some simple algebra to the equations of motion in their final form, suitable 
for explicit analysis: 

6a'' + n{n- 1)13'^ + 6na'/3' = u;e^"+("-^')^ 

^ n(n + 2) 
^ _ (n+p--2)o; 

2(n + 2) ^ ' 

Here primes denote derivatives with respect to r. The specific form of r = r(t) can be 
determined after the solutions are found. 

The equations (|62D are immediately integrable. Indeed, consider the linear combi- 
nation 4a + '^"'^2n+p'^^ f^- using the second order differential equations, it is easy to 



verify that 4a" + ^^^^0^/3" = 0. Hence we can immediately write one first integral 

of m-. 



4„ + M^±iL^^ = C,+C,r (63) 

2n + p 
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where Ci and C2 are integration constants to be determined later. 
To find the other integral of motion, we can define the new variable 



X = 6a + {n- p)P (64) 

The remaining independent second order differential equation becomes 

6n - 4np - n"^ - p"^ ^ 

X = — ue (65) 

n{n + 2) ^ ' 

This is the Liouville equation, corresponding to a particle in 1 dimension moving in 
an exponential potential. This has an easily determined first integral, simply given by 
the conservation of energy: 

X'^ = 2 . ^ ^. uje"" + (66) 
n[n + 2) 

with Eq the energy integral, and where we have defined A = Qn — 4:np — n'^ —p^. At this 
moment, we need to make three observations. First, X is a good independent variable 
only as long as A is nonzero. If it is zero, the variable X degenerates to the previous 
first integral (p3D, up to an overall constant, and so another independent integral should 
be used. This is easy to take into account however, and besides corresponds to a set of 
measure zero in the phase space of solutions, and thus we will not pay it much attention 
here. We will instead focus on the more generic cases where X is independent from 

Second, the curve in the (n,p) plane where A = vanishes is precisely the curve 
defined by ( |59D that separates the traditional potential-free Kasner solutions from the 
more general (but still asymptotic power-law) behavior that we discuss below. In 
particular, the traditional potential-free Kasner solutions apply asymptotically in the 
region A < 0. The full behavior in this region, including the transient regime before 
the power-law dependence on t sets in can easily be discussed by a simple generalization 
of the analysis described below and in Appendix B for the case A > 0. Leaving the 
region A < to the Appendix, we focus on the novel case of A > in the following. 

Third, the integrals of motion C2 and Eq are not independent by virtue of the 
Einstein constraint equation, which is the first equation of ([6^). If we take (^) and 
the first derivative of (|63D, and substitute them into the constraint, we find 

_ ^i^+Pfcl (67) 

4n(n + 2) ^ ^ ' 
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Thus we see that the constant C2 completely controls the dynamics. With this, we 
have essentially reduced the system (|62D to a functional constraint ( |63D and a simple 
1*^* order equation 

X- = 2^^.e- + 3M!cJ (68) 

n{n + 2) 4n{n + 2) ^ ^ ' 

which can be easily integrated when A 7^ 0." 

As shown in Appendix B when A > 0, the exact solutions of these equations yield 
the long-time behavior of the scale factors a{t) and h{t). It is given by 

n(ri+p — 2) 
/ t \ ~ (n+p){2n+p) 

These solutions describe a situation in which the on-brane scale factor shrinks while 
the size of the internal dimensions continues to grow. The exact solutions show that 
for C2 = they are generally valid, while they are asymptotic long-time attractors 
for the generic cases with C2 < (i.e., the case in which both Ha and Hi, > after 
the end of inflation). This means that even if the brane scales initially continue to 
expand, soon after the end of inflation they inevitably start to contract. Note that 
the t-dependence of b(t) is exactly such as to have the potential-dependent RHS of the 
evolution equations scale in the same fashion as the LHS. 

In any case, the main point of this analysis is that the evolution after the initial stage 
of inflation continues into a phase of slow progress of b{t) towards the stabilization 
point, with, under quite general conditions, a simultaneous contraction of our brane 
scale factor a(t). 

The amount of a{t) contraction is controlled by the exponents in (or in the 
case of A < the exponents in (|5^)), and the amount of expansion of b until the 
stabilization point is reached. We can use the asymptotic form of the solutions to 
place an upper bound on the amount of contraction of the brane as a function of the 
evolution of b. We have 



<\j-\ (70) 
a, \ b 



where the parameter ( is given by 



C = !^i^±ZZ^ forA>0 
^ 2{2n + p) 



"if A = a simple modification of this procedure still yields the second integral of motion. 
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c 



3n — ^3n{n + 2) 



for A < 0. 



(71) 
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In these equations recall that A = Qn — Anp — n'^ —p'^, while the effective 4-dimensional 
potential for b varies as V ^ b^^. It is easy to see that the exponents for the a{t) 
and b{t) evolution calculated in the A > and A < regions are continuous on the 
curve A = 0, and thus so is (. Furthermore, note that for a given n the greatest a{t) 
contraction occurs in the A < case. 

3.2 Phenomenology of the era of contraction 

There are a number of interesting consequences of the period of a{t) contraction just 
described. First, during the contraction era, the brane universe looses a number of 
efoldings of the scale factor to the shrinkage. This implies that the early period of 
inflation needs to produce that many efoldings more than the naive minimum needed 
to solve the horizon and flatness problems. Moreover, the contraction of a{t) blue- 
shifts any energy density left on our wall at the end of inflation, and as we will argue 
below, this could, together with the remnant of the de Sitter era Hawking radiation, 
conceivably be the source of reheating of our brane. 

To quantify these remarks, we note that the contraction factor of the brane cannot 
exceed the amount 



where the fact that the final value of b{t) is well-approximated by the stabilizing value, 
bo, has been used. The value of bi is set by the size of the internal dimensions at the end 
of the early period of inflation and at worst is of order the fundamental Planck length 
bi > M~^. (The COBE constraint (|35D typically requires a larger value of bi, and this 
just lessens the amount of contraction.) Thus, the maximum amount of contraction is 
bounded from above by 



The numerical value of this formula ranges between about 900 for n = 2 to about 
2 X 10^ for n = 6. This means that we loose to contraction at most about 7 efoldings 
of inflation for n = 2 up to about 12 for n = Q. Clearly, this is not excessive, and can 
be easily made up for by a slightly longer period of early inflation. 

We now turn to the evolution of the wall and bulk energy densities during the period 
of contraction. First note that the dominant form of energy density at the end of the 




(72) 




(73) 
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inflationary period is the bulk energy of the radion, i.e. the b{t) scale factor kinetic and 
potential energy. The (effective 4- dimensional) kinetic energy density of the radion is 

Pb,KE = nin - 1)M:-^%^{H,Y (74) 

while its 4-dimensional potential energy density is just V{b) as used in the previous 
sections. /^From the exact solutions presented in the previous subsection we find that 
scaling of the energy density of the radion depends on the exponents {k, i) for the 
evolution of the scale factors a ^ t^, and h ~ t^, which in turn depend on whether the 
potential is important or not, in other words the sign of A. The result is 

A > 
A < 



So far we haven't discussed what ends the contraction period. As shown in detail in 
Appendix C, it is a remarkable fact that contraction of a{t) stops and reverses when 
p^aii satisfies^ 

p, = n{n - 1)M:^%^HI = p^au- (76) 

There are two generic possibilities for how this condition may come to be satisfied: 
The first takes the primordial p^aii left over from the inflationary epoch. If there is 
sufficient a(t) contraction, this inflationary p^aii can becomes comparable to pb before b 
reaches the stabilization point, Bq. Then a "Big Bounce" occurs, the contraction stops, 
and a modified expansionary phase that we discuss below begins. This expansion 
finally becomes the usual FRW expansion after b reaches feo- The second possibility is 
that a form of reheating takes place on the wall which is totally unconnected with the 
contraction of a(t), but that again leads to p^aii > Pb- Possibilities in this class include 
the decay of some metastable state on the wall, or the collision of some other brane 
with our brane. 

Consider the first, model-independent, scenario. We are used to thinking of the 
universe at the end of inflation as being very cold, but due to the de Sitter era Hawking 
radiation this is not completely so. What is the initial value of the ratio Pwaii/Pb so 
produced? At the end of inflation, our brane is not empty but inevitably contains 

**The exact analysis we have performed proves that this is the case when the dominant form of 
radion energy is kinetic, and, moreover, as we will show in Section 3.4, necessarily and automatically 
occurs if the radion stabilizes. 



(potential case) (75) 
Pb,PE ~ b~P, Pb^KE ~ b~P 
(simple Kasner) 

Pb,PE--0, 
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a radiation bath left over from the early inflationary de Sitter phase, whose energy 
density at the end of the de Sitter phase has been estimated by Ford p3|] (see also 
g), with the result 

Pu,au,i ^ lQ-\Haj)\ (77) 

depending on the exact nature (conformal or not) of the coupling between the matter 
and metric. Here Haj is the Hubble parameter at the exit from inflation. On the other 
hand, the initial value of the radion energy is of order 

p,, ^ n{n - l)M:^\hr{H,jf ^ QM:+\hjr{Haj)\ (78) 

since inflation terminates when Hi,/ Ha approaches unity. Furthermore due to the slow 
roll conditions, at inflation exit Ha is roughly the same as during the inflationary 
phase. Given these expressions for the initial energy densities, equations (^) and 
(PSP together with the constraint of reproducing the CMB/COBE data lead to the 
interesting relationship 



Pwall 



Pb 



n{n - 1) ^2 f 5_P_ 



(79) 



COBEy 



6 \ p 

Due to the contraction of our scale factor, the wall-localized radiation can become 
important, the ratio pwaiil Pb approaching unity. Using the fact that p^jaii ~ a~^, and 
the relation ([7^) between the evolution of a and 6, we find that it scales as 

P^all ^ _ ^ ^(2n2+4np+p2-4n)/(2n+p) for A > 

P6 ^4 



^-2n+^3n(n+2) 



^V"(2+n)/3 fo^ ^ < 0. (80) 



Given the initial value (|79|) of the ratio we can compute with the aid of the scaling 



laws (pOD how it behaves as the contraction proceeds. Putting in the experimentally 
observed value of 6p/ p ~ 2 x 10^^, and taking a reasonable value of S* ~ 0.01, we find 
that in the simple Kasner A < case, p^aii = Pb before 6 ~ 6o for all n = 2, . . . , 6. On 
the other hand, in the case of the Kasner-like solutions with potential, A > 0, the rate 
of increase of the ratio Pwaii/pb is always slower as a function of b, and in some cases 
{e.g., n = 3, p = 0) the stabilization point is reached before p^aii ~ Pb- We will focus 
on the A < cases for illustration in the following. 

Specifically, using the scaling laws for the evolution of the energy densities, pwaii = Pb 
at a value of 6 = 6/ given by 





-2i/n.(n+2)/3 

(81) 



COBE/ 
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We can now compute the energy density on the wall at the end of contraction: 

/ r . (2(n+l)2-4n^n(n+2)/3) 

P^aiij ^ 7 X 102M,^(6,M,)2"n(n - 1) . (82) 

\ P COBE/ 

In terms of the temperature of radiation at the end of the contraction phase, this 
becomes 



If we use the fact that the typical volume (measured in units of the fundamental Planck 
mass M^) of the internal dimensions during inflation is at most {hiM^)^ ~ 10^, and 
that S6p/p < 2 X 10"'', given COBE and 5" < 0.01, we find that the wall temperature 
at the end of contraction is bounded above by 1.4 x 10~^, 7 x 10""^, 4.5 x 10~^, 3 x 10""^, 
and 0.22 in units of M^,, for n = 2, 6 respectively. 

We compare this to the limits found in [Q], which showed that the temperature of 
radiation on the brane, by stabilization, cannot exceed a certain maximal value in order 
to prevent overproduction of bulk gravitons by evaporation. This temperature, called 
the normalcy temperature T* in 0, was defined as the temperature above which the 
cooling of the wall by graviton production begins to compete with the normal adiabatic 
cooling by expansion. It was estimated to range between about 10 MeV for = 2 to 
about 1 GeV for n = 6. For M^,, ~ 1 TeV, the above expression for the reheating 
temperature of the wall radiation at the end of contraction gives about 1.4 MeV for 
n = 2 up to about 300 MeV for n = 5 and 220 GeV for n = 6. In fact, if we take 
~ 30 TeV for the n = 2 case, as suggested in 0, we can see that the upper bound 
on the predicted wall reheating temperature becomes ~ 45 MeV. Thus we see that 
apart from the n = 6 case we certainly satisfy the normalcy bounds on the reheating 
temperature. 

However, these numbers are only the upper bounds for a number of reasons. First, 
in reality, the evolution is not always Kasner-like. Instead, there are periods immedi- 
ately after the end of infiation when the brane scale factor is increasing, or decreasing 
more slowly than in the Kasner case. As a result of this, the actual blue-shift is less 
than the bounds used above, and this can easily lower the upper bound by an 0(1) 
factor. Second, after the end of contraction, and until the stabilization point b = bo 
is obtained, both a{t) and b(t) expand. Since the normalcy bounds strictly only apply 
after stabilization is reached, and the expansion of a{t) red-shifts the temperature, the 
above estimates of T^aii are actually too high. 




(83) 
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Having found reassurance in these numbers, we can ask if this reheating could be 
not only indifferent, but beneficial for late cosmology. Namely, the bounds we have 
derived above suggest that the reheating of the universe on the wall due to contraction 
may be just sufficient to warm the world enough so that nucleosynthesis can occur 
without hindrance. Indeed, the upper bound on the reheat temperature is always 
above the nucleosynthesis scale, suggesting a very interesting possibility that in this 
picture reheating and particle production could be purely gravitational. In conventional 
cosmological models based on quantum gravity at the scale 10^^ GeV, gravitational 
particle production is typically insufficient for reheating, largely due to the fact that 
inflation occurs early, and that between it and the nucleosynthesis era the universe has 
expanded by many orders of magnitude, diluting the particles produced in the early 
de Sitter phase. In our model, in contrast, the part of the world which contains our 
universe contracts, increasing the energy density and number density on the wall rather 
than diluting them. 

3.3 Bulk graviton production 

We close this section with the consideration of production rates for bulk gravitons 
during the era of radion evolution. There are actually two slightly different issues here 
that we have not so far distinguished: i) the energy density in the (would-be) zero mode 
of the bulk graviton, namely the radion, and, ii) the energy density in Kaluza-Klein 
(KK) excitations of the graviton in the bulk. 

The constraints on the energy density in KK excitations is actually more severe 
than that on the radion energy density, which just comes from overclosure. The reason 
for this is that even though the lifetime of the bulk KK modes is very long, a small 
fraction of them will decay back to photons on our brane, causing distortions in the 
diffuse gamma-ray background [Q] . The diffuse gamma-ray background constraint only 
applies to the excited KK modes, since in general it is only the excited states that can 
decay to dangerous energetic SM states on our brane. (Recall that typically the mass 
of the radion, the would-be zero mode, is > 10^^ eV.) This constraint is more severe 
than the overclosure constraint which applies to both the radion zero mode and all the 
KK excitations together. 

In any case in the next section we will consider the overclosure constraint on the 
total bulk energy, and see that generically there is a problem. Here we want to check 
that the bound on the KK excitations of the graviton is automatically satisfied. (The 
analysis will show in passing that the dominant form of bulk energy will be in the 
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radion motion, and not the KK excitations.) 

First consider the evolution of the projected energy density (of mass dimension 
[Pkk] = 4 as usual) of the bulk KK gravitons in the absence of particle production: 

PKK + ^HaPKK + HbPKK = 0. (84) 

The second term on the LHS corresponds to the usual dilution of massive particles 
(which from this 4-dimensional perspective the KK excitations of the graviton appear 
to be; see for a full discussion). The novel feature of this equation is the third term 
which expresses the fact that as the internal radii increase, the KK masses decrease as 
1/b. Another way of saying this is that the KK mass is really the quantized momentum 
in the internal directions, and this red-shifts as b{t) expands. 

As recognized long-ago by Parker and Zeldovich, and clarified later by many au- 
thors , a time- dependent gravitational field can produce particles from the vacuum, 
by essentially a version of the pair-creation process that takes place in strong electric 
fields for charged particles. Up to order one coefficients that will not be important, the 
magnitude of this particle production per unit time per unit volume is given by H^, 
where H is the typical Hubble constant.^ Thus the equation for the projected KK 
energy density becomes 

PKK + SHaPKK + HbPKK = H\ (85) 

It is easy to solve this equation. Using the substitution pkk = /(^) {^^^ "^^^ ^e 
reduced to 

which, employing the Kasner era power laws, and noting that the resulting power of t 
is less than —1, gives upon integration a result dominated by the early stages 

/ ~ ^ - ^ ^ (87) 

for all power-law solutions. Converted back into Pkk, this says that the final energy 
density of the KK gravitons in the bulk which came from particle production is bounded 
from above by 



Pkkj ^ 



4 / 



a,- \ / 0,- 



^^^'More precisely, as shown in j24j, in the anisotropic case considered here, it is the Hubble constant 
of the contracting dimensions. 
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This differs from our estimate of the final energy density of the blue-shifted wall- 
localized radiation at the end of the epoch of contraction, only in that it is further 
suppressed by a factor of {af /ai){bi/bf), which comes from the fact that the KK gravi- 
tons are red-shifted by the bulk expansion, but only diluted and not blue-shifted by 
the wall contraction. Thus we get 

^ < (89) 

PwallJ Vtti / \bf J 



where in the second line we have used (|70|), and the conservative estimates hi ~ 10M~^ 
and bf = bo. Evaluating this in, for example, the case of the simple Kasner contraction 
with exponents given in ( [7TD leads to Pkk/ Pwaii varying between 3 x 10"^'' for n = 2 
to 1 X 10~^ for n = 6. This shows that the effective temperature of the KK gravitons 
is well below the diffuse gamma-ray bound, even before any dilution necessary to solve 
the radion moduli problem. It also demonstrates that the vast majority of the energy 
in the bulk is in the motion of the zero mode radion pb, rather than in the bulk KK 
modes. This is simply because we have shown above that pf, ~ p^aii is the natural 
circumstance at the end of the epoch of contraction. 



3.4 The era of contraction and stabilization in the Einstein 
frame 

The epoch where the radius grows from its initial small size to its final value can also 
be simply and physically understood in the Einstein frame. (Einstein frame quantities 
will be denoted by an overbar in the following.) Recall that the metric in the Einstein 
frame is related to the one in the string frame via 

9,u = e'^'^g.u, b = (90) 

(throughout this subsection we work in units with M^, = 1). Using this map, it is 
straightforward to relate the quantities in the string frame to their counterparts in the 
Einstein frame. For instance, it is trivial to check that 

Ha = efP (h, - (91) 

Therefore, we can see that even though in the Einstein frame Ha > 0, this can look 
like "contraction" in the string frame, i.e. Ha < 0, provided that /5 is large enough. On 
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the other hand, as long as the radius eventually stabilizes at its final size so (3 — 0, we 
will have Ha > 0. This proves that, in the string frame, there is always a "big bounce" 
as long as the radius eventually stabilizes. 

Moving on to the dynamics, the action in the Einstein frame is 

S^J d'yV^V (^-n+^^^^^^{dW - e-'-~^V(P)^ (92) 
We can parametrize the potential as 

V = b-^W = e-P^J(p) (93) 

We wish to regard / as very slowly varying and will treat it as a constant during 
the radion "coasting". Of course, at the final value of the radion we must have 

m^j'{P)^oj"i(3)>o. 

If we now work with the canonically normalized field 



p = ^n{n + 2)(3 (94) 

the action becomes 

S=j d'y^ (-n + lid^r - e-^^fW)) , X ^ + (95) 

The equations of motion are now just the familiar ones for a FRW universe with a 
scalar field P with an effective potential K//(/?) — e~^^/(/?). Writing the scale factor 
as usual a — e°', we have 



/3 + 3a(3 — xfe 



-x/3 



Suppose we ignore /; then P only has kinetic energy and we trivially find the solutions 

a = ^log(t) 

P = -^log(t) (97) 

Obviously, for sufficiently large x, the effective potential Vgff is falling off so rapidly 
that ignoring / should be a good approximation. We can see what the lower bound 
on X is by substituting the above potential-free solutions into the exact equations of 
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motion; the kinetic energy terms scale like t ^, whereas the potential energy term scales 
hke r^^/^. Therefore, ignoring the potential is a good approximation when 

x>V3^x^-3 = -A>Q (98) 

which agrees with the string frame result. Notice also that this case corresponds to the 
string frame contraction, since 

(Tl ' \ /I Tl \ 

a , ^ = sgn - - < (99) 

2^n(n + 2) J \^ ^3n(n + 2) / 

What about x < \^ or A > 0? In this case, the effective potential does not fall 
steeply enough for the energy density to become kinetic energy dominated. What hap- 
pens instead is that the radion first gets accelerated by the potential till the kinetic 
energy briefly dominates, whereupon it gets diluted again and the cycle repeats. There- 
fore, on average we expect the kinetic and potential energies to be the same. It is easy 
to see that this assumption is self-consistently justified by the equations. Therefore, 
we have 

P = ^ /5 = -log(t) (100) 
Inserting this ansatz back into the equations, we find trivially 

a = ^log(t) (101) 

These solutions correspond to the new, modified Kasner solutions found in the string 
frame analysis. Once again, it is easy to see that these solutions correspond to con- 
traction in the string frame: 

sgn(i/a) = sgn I ^ , ^ ) = sgn(2 - n - p) (102) 

\^ sJn{n-\-2)x) 

while from our solutions in the string frame a contracts when n-\- p — 2 > 0. 

We can now discuss what happens as the radion nears its final minimum (3 = j3^. 
The effective potential Vg// can be approximated as quadratic around this final point 
for {13 — j3^) / 13^ < 1. If the radion approaches this region with kinetic energy sufficiently 
smaller than, or comparable to, the potential energy, then by equipartition it will be 
trapped in the well and will oscillate about the minimum. On the other hand, if the 
kinetic energy is much larger than the potential at the top of the well, (3 will escape 
from the region close to the minimum; its subsequent fate depends on the form of Veff 
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at larger values of /3. It can either escape to infinity, or it may turn back around at 
some large distance away and "slosh" back and forth with very large amplitude about 
It is therefore clear that for x < -\/3, where the kinetic and potential energy stay 
comparable throughout the coasting period, the radion will not overshoot. On the other 
hand, it appears that for x > Vs, where the kinetic energy dominates, the radion could 
significantly overshoot the minimum. This statement has to be qualified, however, since 
all of the above analysis neglects the effect of the radiation energy density left over 
after inflation. Of course, in the Einstein frame, this energy does not appear to be blue- 
shifted, because the Einstein frame scale factor never decreases. However, during the 
kinetic energy dominated era, the kinetic energy redshifts as a^^ whereas the radiation 
energy density only redshifts like a^'^, so the radiation will eventually dominate. Let 
us consider what happens after radiation domination, ignoring the radion potential. 
Since a ~ \/i during radiation domination, the radion kinetic energy ~ a"^ ~ t~^, so 
$ ~ t~^/'^ giving /? ~ t"^/^. In other words, ignoring the radion potential, the radiation 
provides enough friction to stop (3. Therefore, if radiation domination happens before 
(5 approaches the minimum radiation domination prevents overshoot. In all cases 
of interest to us this indeed happens; as demonstrated numerically in subsection 3.2, 
radiation domination takes place for values of the radius h less than the minimum 6o- 



3.5 Stabilization and the moduli problem 

In the previous subsections wc discussed a rather interesting mechanism by which the 
remnant Hawking radiation left over from the de Sitter era could be the source of 
reheating on our brane when combined with the fact that a{t) goes through a period 
of contraction. 

However as we noted above, before this mechanism can be viewed as a realistic way to 
reheat our brane, we must consider the radion field. At the stage where the contraction 
stops and the brane undergoes a Big Bounce, the contribution of the energy density 
of the radion field is still a significant contribution to the total (indeed 0{1)). This 
fact leads to a radion moduh problem. The difficulty is that the radion is so hght 
and weakly coupled that it lives typically much longer than the age of the universe. 
Since its coherent oscillations about the minimum redshift away only as it can 

eventually overdose universe. In order to avoid this, the energy stored in the radion 
must be small, relative to the radiation energy T^, when the universe is reheated to T*: 

Prad* _ Prad* ^ J_ _ PradO ^ J_ ^ 3 X 10~^ GcV fl03~) 
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where we have used p/T^ = const, and Pcrito/TQ ~ 3 x lO^^GeV. Given that T^, is 
bounded between ~ 1 MeV - 1 GeV by normalcy constraints, the energy density in 
the radion must somehow be diluted by ~ 10^^ — 10~^ in order to avoid overclosure. 
Fortunately this implies that, for example, only about 5 or 6 efolds of late inflation 
are needed so solve this problem. In any case this moduli problem is one of late (post- 
stabilization) cosmology, rather than early pre-stabilization cosmology; we will return 
to address it in a future publication. 



4 Conclusions 

We have argued that early inflation when the internal dimensions are still small can 
successfully accomplish all that is required of inflation, including generation of suitable 
6p/ p without the unpleasant introduction of very light or fine-tuned wall fields. Indeed, 
the very fact that the internal dimensions must expand from their initial size close to 
the fundamental Planck length to the stabilization value bo ~ 10~^^^^*'/"cm leads 
to natural inflationary scenarios involving the dynamics of the internal space. The 
resolution of cosmological conundrums such as the horizon, flatness and age problems, 
and the production of the spectrum of nearly-scale invariant Harrison-Zeldovich density 
perturbations with the avoidance of drastic fine tuning of the inflaton mass come 
as a consequence of the evolution of the internal space. Moreover it is remarkable 
that the era of post-inflation brane-contraction that follows this period of inflation 
is harmless, and automatically ends via a "Big Bounce". During the phase of b{t) 
evolution to the stabilization point, the production of bulk gravitons by the time- 
varying metric remains completely suppressed, ensuring that the bulk is very cold 
at, and after, the stabilization of the internal dimensions. The primary remaining 
issue is the radion moduli problem, which is no more severe than in gauge-mediated 
supersymmetry breaking models. Overall, then, early universe cosmology in these 
models is quite interesting! 
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Appendix A: Kinematics of the radion field 



-^ = 6- + 6 - +2n- + nin-l) -] + 6n - + ^— ^, (104) 



With the metric of the form {^) the Ricci scalar is 

2 

a ' \aj ' b ' ^ \b / ' \ab J ' 6^ 

where the internal curvature term is present for n-spheres (k = 1), but vanishes for 
tori (k = 0), and we have ignored a similar curvature term for the large dimensions. 
After integrating over all spatial coordinates we obtain, 

S = J dt{CKE{a, b) - aV(6)), (105) 

and further integrating the d and b terms by parts, the kinetic part of the action for 
the radii, a and b, becomes 

S ^ -Ae" /d*aV (e + - 1) 6„ g)) . (106) 

Note the overall negative sign of these kinetic terms, but also the mixing between d 
and b. 

There is clearly an extremum of the action with d = 6 = 0, when the condition 
da{ci^Vtot{b))\a=ao,b=bo = 0, and similar with da db are met. These imply (for ao 7^ 0) 

Vtot(6o) = 0, and 

VLM = 0. (107) 

This is as one would have naively expected. However, because of the negative sign 
for the kinetic term for the radial degrees of freedom, the stability analysis for such 
static solutions has to be treated with care. The analysis starts by expanding the 
action, Eq. (|106|) , in small fluctuations around the extremum: a(t) = ao + 6a{t), and 
b{t) = 60 + Sb{t). Then to quadratic order, and defining A = 6a/aQ and 6 = 6b/bQ, the 
expansion gives the coupled equations of motion 

(3i,„f-J(f)^(s:0(t)^ « 
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where 

O 71 /f 94-r) / l,\'nO 71/fZ ^'^ 



2 M2+"(6o)" 2 



is the radion mass around the stabihzation point. Searching for oscillating solutions, 
(A, 6) = exp{iflt){Ao, 6o) of the stability equations leads to an eigenvalue problem for 
the frequency fl. Specifically, Q"^ has the eigenvalues = 0, and 

= —^-—uj\ (110) 
n{2 + n) ^ ' 

The zero eigenvalue just corresponds to the fact that ao is a flat direction since, by 
assumption, there is no potential for a. The crucial expression is Eq. ( |110|) , which gives 



the condition for stability of the static solution. In the end, stability just requires that 
the radion (mass)^ be positive as one would expect, and that we can think in terms of 
a total potential V{h) that one can minimize to find the stable static solutions for the 
size of the internal dimensions. 

The equations of motion for a{t) and h{t) derived from the action Eqs. (|105|) and (|106|) 



are, after some algebra, those given in the text, namely (0), (|TT]) and the constraint 
(0) which comes about from the well-known property that the total energy in GR is 
zero. (This constraint can be derived by carefully working in terms of the lapse and 
shift functions of the canonical formalism.) If matter on our wall is also included then 
they become those given in (p^): 

V + p 



QHt + n{n - l)Ht + QnHaHi, 



^ / N ..2 I ( 2V b dV p-3p 



b ^ ' " M,"+26" \n + 2 n{n + 2) db 2{n + 2) J 

'i.OH^^^HH _ ^ ( h dV n-2 p + (n - l)p \ 

a " ^ " M::+%^\2{n + 2)db 2{n + 2) 2(n + 2) / 

p + 3Ha{p + p) = Q. (Ill) 

Note that in these equations the effect of wall-localized matter is just some extra 
contribution to the a{t) and b{t) scale-factor evolution. Of course the energy density 
on the brane in general distorts the geometry of the internal space (as does that present 
on other branes that may exist in the bulk), but as far as the overall properties and 
evolution of the zero mode size modulus b{t) of the internal space is concerned, it is 
correct to treat the energy density on the wall as just averaged over the whole space, 
as done on the RHS of these equations. 
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It is useful to summarize some basic properties of the evolution equations. Since the 
total potential energy U{a, b, ip) is given in terms of the effective 4-dimensional potential 
energy density V{b,ilj) hj U = a^V{b,tp), a uniform bulk cosmological constant is 
represented hj V = 6" A. Substituting this form into the RHS of these equations we 
see that a positive bulk cosmological constant term has the effect of wanting to increase 
both a and b if the evolution of a, b are studied close to zero. This is not inconsistent 
with the stability criteria derived above since this was explicitly the stability analysis 
around a stationary point of the equations with non-zero values of both a and b. 
Indeed, the stability analysis can be derived directly from the equations of motion, as 
it must. Concretely, if we expand around a point (a*, 6*) with V{b^) = V'{b^:) = 0, 
then linearizing the equations of motion gives 




V" 

n{n+2)Mf^ ^ 




(112) 



which exactly reproduces the previous stability analysis, in particular the requirement 
V"{b,) > for a„b, > 0. 

Another basic property that is useful to keep in mind is the effect of some small 
amount of wall-localized matter on the position and mass of the radion. Linearizing 
the equations ( p7[ ) with such matter in the shift 6b = Sb/bo around the stabilization 
point we find 

A^-^Aff 1 (blV" n{p-3p) \ {p-3p) 

Sb + SSbH^ - -(^ + 2)6W+2 + 2 ) + 2(n + 2)6W+^' ^ ^ 

This shows that wall- localized matter has two effects o the radion. First, the 6b- 
independent term on the RHS shifts the radion from the stabilization point by a small 
amount, and second the (p, p)-dependent in the parentheses on the RHS shifts the 
eigenfrequency of oscillations of the radion (or equivalently the radion mass). Both 
these effects are to be expected, and are harmless for wall-localized matter densities 
and pressures p,p ^ M^. 



Appendix B: Kasner-like solutions with potential 

The solutions of (0) take different form, controlled by whether C2 vanishes or not 
and the sign of A. Let us first consider the case A > 0. Then, when C2 = 0, the 
equation (^) simplifies to X' = ±^^2^^^^^^; exp(X/2), which can be integrated to 
give exp(X) = 72 • Here we have removed an additional integration constant 
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by a time translation. It is straightforward now to use this and (^) to determine the 
solutions for a and (3. Since a = (2n+PKn-p)Ci _ n(n+p-2) ^ ^ ^ _ 3(2n+p)Ci ^ 

we find 



n(n+p — 2) 

(2n + — p)Cl \ / ACJ \ n(n+p-2) 



^ = ^"Pl 4A AM^' " ^ 



2n+p 
A 



r A 



We can now use (|6ll) to transform (|114|) back to using the comoving time on the brane. 
Integrating (|6lD , we find 



r ~ t (i+p){2n+p) _ (115) 

Note the important feature of this map, that it maps the comoving future t ^ 00 to the 
origin of "time" r and vice-versa. Thus to discover the correct long-time asymptotics 
of the solutions we must extract the small-r behavior. This is true for the general case 
C2 7^ as well. 

In any case, for C2 = 0, and after appropriate rescalings, we find the behavior of the 
scale factors a{t) and b{t) to be 

n(n4-p — 2) 
t \ ^ (n+p){2n+p) 



These solutions describe the cases where the brane length scales shrink while the radion 
continues to grow. This means that even if the brane scale factor a{t) continues to 
expand by inertia after the end of inflation, soon after a{t) starts to contract. Moreover, 
the behavior (|116|) is an asymptotic attractor for the generic cases with C2 < and 
A > 0, as we will now discuss. 



Let us now consider the case when (72 < 0. Again, note from (|63D that this corre- 
sponds to the situation where initially both Ha and Hi, are positive, when inflation has 
just ended. The equation ( |68D can be solved in this case by the substitution 

which reduces to the differential equation 'd' = ( i6n"n+2) ) ^ ^2, whose solution is ^9 = 
( i6n(n"+2) ) ^ ^2T, after the appropriate choice of the coordinate origin of r. Therefore, 
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the solution of (|68|) is 



3(2n + p)2 

3(2n + p)^Cl I6n{n + 2) 



sinh 



:il8) 



Note that the RHS is a positive semidefinite function, guaranteeing the reahty of the 
metric, as required. Also note that the solution (|118| ) is defined in open interval 
r G (0^,00) and (— oo,0~), by time-reversal. Using this equation and (|6^), we can 
extract the solutions for a and (3: 



8Au 



3{2n + v?Cl 
3(2n + p)2C| 



n(n+p — 2) 
A 



e 4A (C.l+G2Tj gj^J^ s 



3(2n + p)2 
^ 16n(n + 2) 



C2r\) 



8Au 



^exp(-^(Ci + C2r)) 



sinh' 



, 2n+p 



3(2n+p)^ 
16n(?i+2) 



|C2r|) 



;ii9) 



where the argument of the hyperbolic sine is taken to be positive, to make sure that the 
solutions remain real. In general, it is not possible to explicitly determine the integrated 
form of the gauge transformation (|6TD . However, the asymptotic limits r — 0, 00 are 
easy to deduce, and are sufficient for our purpose here. 

Taking the r — limits of (|119|) we find that (|119|) reduce precisely to ( |1 14|) , 
implying that the integrated form of (|6TD approaches (|115|) . Hence in the case C2 < 
the scale factors again approach (69) as t —>■ 00, which are therefore the appropriate 
future attractors in all cases A > 0, C2 < 0. 

It is also amusing to consider the short t-time behavior of the exact solutions in this 
case. We find that as r -h> cxo. 



exp(^ 
exp(^ 



{2n + p)C2T 
4A ' 
{2n + p)C2T 
2A 



n — p — n{n + p — 2 
(3 - (2n + p) 



n{n + 2)' 



nin + 2) 



Using this, we can integrate (|6T|) and find 

■{2n + p)C2T 



exp 



(■ 



{n{n — p + Q)\ 



3(n + p 



(120) 



(121) 



4A ' n(n + 2) 

and upon substituting this back into (|120|) , we can show that the solutions approach 
from below the asymptotic expressions 



J. 3+Y^3n(n + 2 
3(n + 3) 



39 



The dependence on the parameter p apparent in { \12(\ ) has completely disappeared. 
This can be most easily seen by taking the derivatives with respect to p of the powers 



which are obtained by substituting (|121|) into ( 120 ), noting that they are identically 
zero, and then setting p = to obtain ( |122| ). Also note that the resulting powers are 
identical to those found in the simple Kasner case (0), with lower sign taken. Thus 
in the very early time limit, this suggests that the universe was expanding while the 
radion was decreasing. However, this phase is cut out by a stage of inflation, and in 
fact only a very short portion is retained where both a and b are growing for short 
time, to match onto the post- inflationary era. 

Now we consider the last case, A < 0, which we argued in the text should be 
described by the traditional potential-free Kasner solutions in the long-time limit. We 
can see this explicitly by examining the exact solutions once again. First note that, 
in this case, C2 must be nonzero, as can be immediately seen from (|68|) . The exact 
solutions can be found by the substitution 



-X 



8\A\u 



■ cosh^(^?) 



(123) 



3{2n + pYCI 

which is analogous to (|117|) . The equation for -d is the same as before, and hence the 
solution is 

SlAlcj 



-X 



cosh 



?>{2n + py 
^ lQn{n + 2' 



■C2T 



(124) 



3(2n + p)2C| 

The similarity of this solution to (|118|) allows us to extract the expressions for a and 
/? quite easily. We find 



SlAlcj 



3(2n + p)2C| 
?,{2n + pfCl 



^ I /I '- (2n+p)(n-p) 



{C1+C2T) 



cosh' 



ri(ri+p — 2) 



3(2n + p)2 
^ 16n(n + 2) 



\C2r\ 



SlAlcj 



2n+p 
|A 



exp(^- 



3(2n+p) 
2|A| 



(Cl + C^t) 



2 2n-\-p 

cosh 1^1 



3(2n+p)2 
16n(n+2) 



\C2r\ 



(125) 



Note however the important difference between these solutions and (|119|) . In this case, 
the solutions are defined on the whole interval (—00,00). By considering the limits 
T —>■ ±00, we can see that the solutions (|125| ) in fact interpolate between the simple 
Kasner solutions (|52D, where for C2 < it starts out with powers given by ( ^4]) with 
lower sign and transmutes due to the intermediate potential-dominated region into (^2]) 
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with powers (|5^) with upper sign taken. In reahty, the far past of these solutions is 
cut out by the inflationary era, and again only a very short portion where both a and 
h are growing is retained, which is matched onto the post-inflationary era. All of these 
solutions therefore flow towards simple Kasner solutions with the powers with upper 
sign in (|5^), which are the appropriate future attractors. 



Appendix C: Exact Solutions for the Big Bounce 

To see that bounce behavior does indeed occur in the presence of radiation on the wall 
we go back to our equations of motion but now set the potential to zero and just keep 
the wall- localized radiation terms on the RHS: 

r2 , 1 Mt2 , ^„ TT TT P 



mi + n{n - \)Ei + ^nHaBi, 



6 ' ' " " " 2(n + 2)M,"+26'^ 
- + 2H^ + nH,H^ - ^^^^^^ + 2) 

p + 3ff,(p + p) = 0. (126) 

These can again be solved exactly, and provide a good approximation for the exit from 
the Kasner-like phase of a(t)-contraction. 

Set p = pI?> for radiation. The solution for the radiation energy density is as usual 

5M"+2 , ^ 

p = (127) 

(with some normalization B) and define variables a and (3 by 

a = aoe° 

b = boe^ (128) 



where oq and bo are constants, determined by the values of a and b at the end epoch 
when the radion potential becomes small compared to the wall-localized radiation. 
Using this, the equations of motion can be rewritten as 

6a^ + n{n - 1)$^ + 6na$ = Ae~''^-^'' 
fi + p{3a + np) = 

ci + d(3d + n/3) = 4e"^"""^ (129) 
6 
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where A = B/ (oo^o). If we change the time variable, defining the new time r as in eq. 

my- 

dr = -e-3°-'^^dt (130) 

where again the specific form of r = r(t) can be found later when solutions are deter- 
mined, we can rewrite the equations of motion as 

6a'2 + n{n - 1)(3'^ + 6na'(3' = ^6^°+"^ 
(3" = 

a" = -e^^+'-P (131) 
6 

where the primes denote derivatives with respect to r. The (5 equation immediately 
gives 

(3 = Ci + C2T (132) 

Hence note that if if;, > initially then the constant C2 < 0. Also note that given our 
conventions we again have future t infinity map to r = 0. 

For the purpose of examining the solutions to these equations we distinguish two 
cases: 

(i) C2 = 0; like before, this solution turns out to be a future attmctor of all solutions] 

This case has a very simple analysis: C2 = implies /5' = but also /? = 0. Hence 
we can forget about the r coordinates and immediately work in our original brane-time 
t. So 6 = means b = bo = const and so the constraint equation gives 

= — 7 (133) 

which is immediately solved to give us the radiation dominated universe, a ~ Vt. This 
is precisely as we would expect from an analysis in the Einstein frame. 

(ii) C2 7^ 0; this case corresponds to more generic solutions. In particular we will 
see that we can have initial conditions H}, > 0, Ha < 0, which evolve to Hf, = 0, 
Ha > 0, in other words a "bounce" solution. 

To solve the equations in this case, first define 

X = 2a + np, (134) 
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in terms of which, the second order differential equations produce 

X" = (135) 
This is again the Liouville equation, with the first integral (conservation of energy) 

X" = + Zo (136) 

where Zq is another integration constant. If we eliminate from the constraint 
equation by using the last equation, and then simplify the result by using the second 
equation, we find that the constraint reduces to 

Z„ = !!(!!±llc| > (137) 

Thus the constant C2 controls the dynamics, and the whole system has collapsed down 
to two simple equations, one of which is already solved: 

2A ^ nin + 2) , , 

X'2 = — + ^ g ' Cl (138) 



Note that the solution to the X equation is, 

^ _ n(n + 2)Cl 1 



(139) 



12 

as may be checked by substitution. Using the formula for X = 2a + n(3 and 13 = 
Ci + C2T, we find that 

V ; sinh(V^IC.Ir)' ^ 

where since C2 < in the cases of interest we have taken the appropriate branch of 
the square-root such that a = aoe" is positive as it must be. 

This expression is already sufficient to show that we get a bounce behavior for a{t). 
Recognizing that t = corresponds to r ^ 00 while t —>■ 00 corresponds to r ^ 0, 
simply plotting as given by ( |14(J| ) shows that a(t) goes through a bounce. 

It is instructive to see this in detail. First, consider the limit r — 0. Clearly, 
P = Ci + C2T ^ Ci = const. Also, 

/Se-nCiX 1/2 ^ 

— - (141) 
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Thus exp(3a + nf3) — r ^, and so we have 



dr 

dt (142) 



r3 



or, t ~ r ^ which gives us r ~ 1/vt, and t ^ oo maps to r 0^ as claimed. Back in 
the formula for a, this gives 

a ~ exp(Q;) ~ — ~ (143) 
r 

and this is precisely a radiation-dominated universe at late times! So indeed, the 
solutions with C2 = are late time attractors. Some further analysis shows that in 
this limit, (and in the approximation of ignoring the stabilizing potential) b will tend 
to a constant logarithmically. The end result of the analysis is that for large t these 
solutions show that a is expanding asymptotically as ^/t. 

Let's now look at the other limit, r 00. Also, recall that C2 must be less than 
zero: C2 < 0. We have 

b ~ exp(/3) ~ exp(-|C2|r) (144) 



and 

a ~ exp(aj ~ exp I — - — 
These give the relationship between r and t 



(n- y'n(n + 2)/3) ) . (145) 



dt expl^^in- J3n{n + 2))\dT (146) 



Now note that {n — ^J3n{n + 2)) < for all 2 < n < 6, and therefore as r — +00 we 
have t — > 0+ as 

'-\C2\r, 



t^expl^ {^3n{n + 2)-n)j (147) 

as claimed. Finally, in the expression ( |145| ) for a(t), note that (n — ^n{n + 2)/3) > 0, 
so as r ^ 00, the scale factor a{t) is again going to infinity. Re-expressed in terms of 
t we see that as t initially increases a{t) is initially decreasing. 

If one carefully considers this limit the power law behavior of a and b in terms of the 
time t precisely corresponds to the late-time Kasner solutions we found in Appendix B. 
Most importantly, we see that as r decreases from cxd towards 0"*", the radion b is 
initially increasing and a is initially decreasing, the initial conditions that we require. 
But as we have seen from the analysis in the r — limit, this goes over to b growing 
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logarithmically and a increasing as 1/r. This implies there must indeed have been a 
Big Bounce in between! 

Finally, the bounce occurs when d = or equivalently Ha = 0. This happens when 



n(n - l)Hl = .^""f. ■ fl48) 



Hence, when the bounce occurs the radion kinetic energy is comparable to the wall 
energy density. By continuity, it should be clear that the generic qualitative features 
of these properties would remain true even in the presence of stabilizing potentials. 
The main conclusion of this analysis is that the Big Bounce is the future asymptotic 
attractor of all post inflationary solutions with wall radiation, and hence the exit from 
the contraction on the wall will occur naturally. 
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